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We report on lattice results of the form factors for semi-leptonic decays of the Dj meson to Tj and 
T]', with n/ = 2+ 1 configurations. The calculation contains disconnected fermion loop diagrams, 
which are challenging to calculate on the lattice. Our result shows that the disconnected diagrams 
give significant contributions to the form factors. 
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1. Introduction 

About 15% of the decays of the meson are semi-leptonic. There are by far dominated by 
decays with rj, rj' and <p mesons in the final state. The decays into rj and rj' contain interesting 
flavor physics, in particular, the mixing of rj and rj'. Investigating the corresponding form factors 
helps to understand the magnitude of the gluonic contribution to rj' (see for example). The 
mixing can be studied through both D and semi-leptonic decays. However, the latter is easier to 
access in experiment, since it is a Cabibbo allowed process. 

The Ds — )• T] ( ' semi-leptonic decays are characterized by the following two form factors, 

/o(^2)and/+(^2). 



{ll^'\k)\V^{q^)\D,{p))=f+{q') 



+ Mq') ' , ' q', (1.1) 

q 



where is a vector current and Mo, and M^(i) are the masses of and rj , respectively. Together 
with the CKM matrix, one can relate the form factors to the decay width and thus the branching 
ratio. Experimental results for these decay modes are starting to appear, see [Q, ||]. The matrix 
element (l.h.s) is what we can obtain from the lattice. We will focus only on the scalar form factor 
fo{q^), which can be related to the scalar matrix element [Q] 

where S = sc, and nis and nic are the s- and c-quark masses, respectively. This relation has an 
advantage over eq. ( |l.l[ ), in the sense that it does not require the current renormalization. 

Currently, only a prediction of the form factors from the light cone sum rules is available 
A first principles calculation on the lattice is desirable. It also provides an interesting playground 
for quantum field theory, since physics relevant to rj' should contain non-trivial effects from the 
chiral anomaly as characterized by the Veneziano-Witten formula 

However, a lattice calculation of the decay form factor is technically challenging, as it re- 
quires the evaluation of disconnected fermion loop diagrams. In order to obtain the relevant matrix 
elements, we need to compute the following three point functions shown pictorially: 




where the solid lines represent fermion propagators. The first and the second terms are the con- 
nected and disconnected fermion loop diagrams, respectively. Note that the disconnected one is in 
fact connected by gluons. We expect to see the effect of the anomaly in the pseudoscalar flavor 
singlet sector and this is only possible with the disconnected terms. 

To obtain a fermion propagator, we need to invert the Dirac operator, a matrix of typically 
O(IO') X O(IO'). Many of the quark propagators in the three point function appearing in eq. (^3|) 
only require a few columns of the inverted matrices. However, the disconnected diagram requires 
the whole inverted matrix (i.e. an all-to-all propagator) for the light quarks (/) and thus is compu- 
tationally much more expensive to calculate. Note that the disconnected contributions are summed 
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over the three light flavors, which enhances the contributions by roughly a factor of three. This is 
another reason to expect that they may be large. 

Although the calculation is challenging, it is still feasible |j8|]. We use stochastic estimations of 
all-to-all propagators in the disconnected fermion diagrams. In addition, we also use the stochas- 
tic method for the charm quark propagator in the connected diagram (denoted by a blue line in 



eq. ( |1.3| )) [^. This allows us to access many different momentum combinations, which helps to 
improve the signal by averaging over rotationally equivalent momenta. The stochastic estimations 
introduce stochastic noise in addition to gauge noise, and we combine several noise reduction tech- 



niques, such as the standard low mode averaging and the truncated solver method (TSM) [10] 



To carry out the lattice calculation, we use QCDSF rif = 2 + \ configurations [ ]11[ , |12| ]. The 



strategy of these sets of configurations is to keep the flavor singlet mass average of the three quarks, 
^{mu + md + ms), constant. Starting from the flavor SU(3) symmetric point, mi,,d = rris, rriu^ is 
reduced and lUs is increased. This setup is ideal to study flavor physics in the SU(3) flavor basis. 
So far, we have analyzed 939 configurations at the flavor symmetric point, corresponding to pion 
and Kaon masses Mj^ = Mk = 450 MeV, and 239 configurations with Mjc = 348 MeV and Mk = 
483 MeV. The physical spacial extent of the lattice is roughly 1.9fm, and the lattice spacing is 
a ~ 0.08 fm for both sets of configurations. The stout link non-pert urbatively improved clover 



action (SLiNC) [|13|], for which the discretization effects are removed to 0{a), is used for both the 
dynamical u, d and 5-quarks and the partially quenched c-quark. 

In the next section, we extract the masses of rj and rj', and a mixing angle between the SU(3) 
flavor singlet-octet basis and the physical states. We then present our result for the scalar form 
factor /() in section |3[ 

2. Extracting rj and rj' 

Since the configurations are well suited to study the effects of flavor symmetry breaking, we 
start with the following SU(3) basis to describe Tj and rj': 

rjs = —^{uii + dd — 2ss), rji = —^{uu + dd + ss). (2.1) 

V 6 V 3 

The physical rj and tj' should be a mixing of the above octet state (rjg) and the singlet state (t]i): 
77 =cos0T]8-sin0T]i, T]' = sin0T]8 + cos0T]i. (2.2) 



We use a single angle parameterization and do not include a gluonic state in eq. (2.2), since we 
expect the gluonic contributions to be small. In principle we can resolve the contribution but is 
beyond the current analysis. The interpolating operator ff^^,^ that will be used creates an overlap 

with the T](') state sufficient to obtain the target matrix elements. 

To extract the rj and rj' states from the lattice data, we need to diagonalize the following 2x2 
two point correlation matrix: 
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where and are interpolating operators for rj^ and rji, respectively, and t is the time separation 
between source and sink. The operators at the sink are projected to momentum k. Here we use 
smeared interpolating operators (Wuppertal smearing [ 14 1 ) to reduce the effects from the excited 



states. The diagonalized basis gives the operator 
point functions give the spectra for ?/a » 1 : 

17 ^|2 



for physical T] , and the diagonalized two 



c; 



^r\,k' 



2E. 



IZ 



IE. 



exp(- 



(2.4) 



where Z, 



(t]('^(^)|^^(,j(^)|0) is an overlap factor between the r\^'^ state with momentum k 



and a state created by the interpolating operator. E ,i) ^ is the energy of the state. Setting k 



0, 



we obtain the masses. Note that each element in eq. (2.3) contains both connected and discon- 



nected fermion loops. For example, the T]8 
diagrams: 



two point function contains the following fermion 



i{t;k)ffl{0)) 



I 

I I 



s 
s 



-2(#C3 Q#)-2(#C3 C=#)+2(#C3 Q#)+2(#C3 C=# ) 

where I stands for light (i.e., u- and d-) quarks. In figure |l|, we show the effective mass plot 



aMeff ( t + 



l0£ 



(2.5) 



for T] and r\' together with that of n. The figure demonstrates our ability to separate the T] and r]' 
states. 

The diagonalization presented so far gives a mixing of the interpolating operators &\ and i^g, 
which depends on the smearing of our choice. To obtain the physical mixing, we replace the 
smeared operators at t in eq. (2.3) by local operators. The smeared operators are still used at f = 
so we need to use Z^(i) from smeared(source)-smeared(sink) two point functions to normalize the 



matrix elements. Setting the momentum k 



0, we can extract the following matrix elements: 



Modulo renormalization factors, they are proportional to the decay constants in the SU(3) flavor 
basis: for example, Zga^ = M^f^, where Zg is the renormalization factor and is the decay 
constant of T] through the octet. We use the following ratio to obtain the mixing angle d, 



^Q|^local| 



(0|<=^'|t]') 



(2.6) 



tan^e 



(2.7) 



(0|^^'=^'|77)(0|^|°'=^'|77')' 

for which the renormalization factors cancel. 

Figure ^ shows our preliminary values for the mass (left panel) and mixing angle (right panel) 
together with the experimental values and other lattice results. We findM,j = 513(ll)MeV, My^i = 
750(130)MeV, and 6 = -8.3°(2.8) at = 348MeV. Note that our resuh has no mixing by 
definition at the SU(3) flavor symmetric point: T] = Tjg and Tj' = Tji. Currently, we have not yet 
included disconnected two point functions for this ensemble, which means M-q = M-^i 
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Figure 1: Effective mass plot for n, rj, and rj' for the Mj^ = 348 MeV ensemble. 
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Figure 2: Summary of lattice results for the masses of rj and rj' (left panel) and the mixing angle (right 
panel). Included are result from the BRC/UKQCD [|l|], UKQCD [|l6^ HSC and ETMC ^ col- 
laborations. For the mixing angle, t'wo experimental values by KLOE [|l9|] (•, o) are from radiative decays 
(p, o, 0) — >^ {n^,ri,ri')Y, rj' — > py and(7r'^, J]') 77 ; they assume zero (•) and non-zero (o) gluonic contri- 
butions to rj' state, respectively. The other experimental value (A) is from BES [pO|], ■which uses charmonium 
decays only. Data ■with * use both local and fuzzed (smeared) operators and ** fuzzed (smeared) operators 
only. 



3. Decay Form Factors 

Above we have constructed suitable interpolating operators for Tj and rj', so it is straightfor- 
ward to obtain the matrix element we need from the following three point function: 



cf'-'''\o = (^(')(^/;W;^)<(0) 

-J" -Tzr = exp 



(3.1) 



Here, p = k + q and i^d, is an interpolating operator for Dg. We use a similar parameterization 
to eq. {2A) for the two point function {&dX^'^p)^d (0))- combining eqs. (lA), (3J) and 
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(^D,(^; p)^D,(0)), we obtain the matrix element {ri'-''>){k)\S{q)\D,{p)) . 

In figure ^ we plot our preliminary results for fo{q^) for the SU(3) flavor symmetric case. A 
fit to the functional form foiq^) = is also shown. We find /o(0) = 0.75(3) for D, fv/T] and 
/o(0) = 0.52(5) for D, /"v/Tj', at Mj, = 450 MeV. 

To calculate the three point functions, we included all disconnected fermion loops. As men- 
tioned in the previous section, we have used the approximation Mjj = M^,/ to compute the two point 
functions, which allows us to use different 4ep = tf — ti in eq. (3J) for the connected fermion loop 
contribution (4ep = 24fl:) and the disconnected one (fjep = 8«)- A larger 4ep suppresses possible 
pollution from excited states, but leads to larger statistical errors, which is particularly problem- 
atic for the disconnected diagrams. Using different fjgp for the connected and disconnected three 
point functions circumvents these problems in this case. A full analysis with Mrj t^M^,/, which 
could affect the Dg — >• Tv/T]! form factor, is currently under investigation. Due to flavor symme- 
try, the Ds — > /V/Tjg form factor is identical to D — )- IviK and D —?- IviK at the symmetric point 
{M^ = 450MeV). 

Interestingly, we find that the contributions from the disconnected fermion loops are signif- 
icant, see figure ^ where both connected and disconnected contributions to C^'^^ (t) are plotted 
separately. The figure also indicates that the main source of the statistical error comes from the 
disconnected contribution. 

The results presented here are still preliminary. In particular, the effects of excited states of 
the Ds and possibly T], T]' are still under investigation. The effective mass plot (figure suggests 
that the excited state may contributes up to f/a < 5. Since our fit region for the three point function 
is at most 1 < t/a <1 for the disconnected diagrams, there may be some pollution from excited 
states. Note that the effects of an excited state (e.g. rj*) to the two point functions is suppressed 
by a square of the overlapping factor. For example, including the leading order of the excited 



contributions eq. (2.4) becomes 



c: 



,diag,?7 



t,k) 



\Z -|2 



2E, 



ri,k 



exp(-£^/) 



ri*,k' 



^r\,k 



IZ E - 



exp 



-(E ,7-E r)t 

^ r\* ,k r\,k' 



+ ••• 



(3.2) 



The term |Z^*/Z^p appears as a correction. On the other hand, the pollution to the three point 
function is suppressed by only a linear factor: 



C 



(0 



4E„rEn 3 



exp 



7 - F - 

1 H — — exp 

7 F ^ 
^r],k ^ri\k 



ED,rp(t-ti)-E^^r,j^{tf-t)\ {ri(k)\smD,{p)) 

{71*ik)\S{q)mp)) 



+ ^|^exp [-(Eo^^^p-Eo^rp)it-U)] 
Zd„p Ed*,p 



{7i{k)\smD,{p)) 

{r]{k)\smDt{p)) 
{r](k)\S{q)mp)) 



+ ■ 



(3.3) 



where the correction starts with the Zri*/Zri term. It is important to remove these effects or to 
estimate their magnitude, which may be one of the main sources of the systematic error in our 
calculation. 
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Figure 3: Preliminary results for the scalar form factor foiq^) for the SU(3) flavor symmetric ensemble 
(M:^ = 450MeV). Also included ai-e results from light cone QCD sum rules (LCSR) [||. 



p=000, q=000, k=000 



O 



5e-07 
4e-07 
3e-07 
2e-07 
1e-07 


-1e-07 
-2e-07 
-3e-07 



full 

connected only 
disconn. only 



4 5 

t/a 



Figure 4: Connected and disconnected contributions to the three point function C3 (t) with p — q = k — 
and M,i = 348 Me V. The meson is created at f /a = and the Tj ' is annihilated at f /a = 8 . 



4. Conclusions 

We calculated the semi-leptonic decay form factors for — )■ / V/T] and — /V/T]' using 
QCDSF Hf = 2+1 lattice configurations. This is the first lattice result of decay form factors 
which includes the fermion disconnected loop contributions. It turned out that the disconnected 
fermion loops give significant contributions to the form factor. The result is still preliminary but 
promising: the scalar form factor at zero momentum transfer /o(0) can be obtained with 10-15 % 
precision. 

We are planning to use larger lattices and lower quark masses. Another target in addition to 
fo{q^) is the other form factor f+{q^). We are also interested in calculating the — )• Ivi^ decay 
form factor that also contains a disconnected quark line contribution and the same calculation 
techniques are applicable. 
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